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The idea of saturation of parton densities in small x physics is briefly 
| introduced. Some aspects of saturation are described, mainly focusing on 

QO . the status of our knowledge on the non-linear equations describing the high 

parton density regime. Implications of saturation ideas on the description 
of nuclear collisions at the Relativistic Heavy Ion Collider are discussed. 

q ■ PACS numbers: 12.38.-t, 24.85. +p, 25.75.-q 

1. Introduction: small x and saturation 

The BFKL equation resums gluon ladders taking into account all 
r~| , leading contributions [a s In (l/x)] n (LLl/x), with a s = a s N c /7r and x the 

fraction of momentum of the hadron carried by the parton. Although it was 
originally an attempt to compute the high-energy asymptotics of QCD and 
rS [ to justify Regge Theory, it turned out to be an evolution equation in 1/x for 

the so-called unintegrated parton distributions (used in /cr-factorization [2] 
to compute inclusive particle production at scale Q, Aqcd <C Q <C E cm ), 
which behave oc x _Qs41n2 . Experimentally |3j xG{x,Q 2 ) oc x~ x , with A ~ 
0.3 sizeably smaller than predicted by BFKL for reasonable values of a s . 

Both /cr-factorization and the BFKL equation are valid in a low parton 
density, linear regime. At small x (equivalent to large energies for a fixed 
scale) parton densities become high. Then the idea of saturation of parton 
densities (see [3]) becomes unavoidable, as parton fusion balances parton 
splitting if A^ u (xl/g. It can be alternatively formulated in terms of the S- 
matrix at fixed impact parameter b becoming black, \S(b)\ = 0. Saturation 
constitutes a new regime of QCD, in which old ideas (pomeron interac- 
tion, multiple scattering, coherence arguments,. . .) are expressed in a new, 
QCD language. It also offers a link between small x physics and heavy ion 
collisions: in nuclear collisions at high energies large partonic densities are 
expected, both due to the high energy and to the ^4 1//3 -enhancement coming 
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from the nuclear size. The understanding of the initial state in a nuclear 
collision is crucial to get a coherent picture of the eventual thermalization 
of the system and creation of Quark Gluon Plasma. 

Our present understanding of the different regimes of QCD is summa- 
rized in Fig. ^ At low Q 2 we have the confinement region. At large Q 2 and 
not too high 1/x we are in the low-density region, where the usual evolution 
equations can be used: BFKL in In (1/x) and DGLAP in InQ 2 . For both 
large In (1/x) and InQ 2 we are in the DLL regime, where a first non-linear 
correction has been proposed, the GLRMQ equations [3]. Finally in the 
high-density regime, separated from the low-density one by a line which 
defines the saturation scale Q s (x), a non-linear generalization of BFKL is 
the Balitsky-Kovchegov (BK) equation [EJIIj- 




Fig. 1. Regimes of QCD in the In (l/x)-\nQ plane. 

I will present in the next Section the phenomenological models imple- 
menting saturation and our current theoretical understanding of this new 
regime. In Section 3 I will analyze properties of the solutions of the BK equa- 
tions, and in Section 4 the application of saturation ideas to high-energy 
nuclear collisions, presently studied at the Relativistic Heavy Ion Collider 
(RHIC) at BNL. Finally I will draw some conclusions. Due to space limita- 
tions, many interesting topics like NLL developments 1 , non-linear equations 
in the collinear approach, relations with other realizations of collectivity 
or many examples of observables analyzed, will not be discussed (see e.g. 
[U IH1 E] and references therein) . 

1 NLL effects are expected to be important at iarger rapidities, Y > l/a 5 J s , than 
unitarity corrections, important for Y > In (l/a s )/a s . 
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2. Present realizations of saturation 

2.1. Phenomenological models 

The most compelling indirect evidence of saturation comes from the 
phenomenological success of some models containing saturation ideas when 
confronted to HERA data. The most commonly used is the GBW model 
jlUj . It is formulated within the color dipole model, in which the interaction 
of the virtual photon with the hadron or nucleus is described as a convolution 
of the probability that the photon fluctuates into a qq pair of fixed transverse 
size r with the dipole-hadron cross section 2 . For the latter, the GBW model 
provides an ansatz for the scattering amplitude N: 




which implements the unitarity limit, a q ^~ h (r,b) = 2 Re[l — S(r,b)] < 2, 
in a very simple way. With oq ~ 20 mb, xq ~ 3 • 10~ 4 and A ~ 0.3, this 
model gives a reasonable description of all HERA data on F2 of the proton 
for x < 0.01 (Q 2 < 450 GeV 2 ). It has been extended to included DGLAP 
evolution and to describe diffraction at HERA and is widely used for 
phenomenology [T^J. It implies a [r = Q 2 /Q 2 (x)]-scaling observed 
at HERA and also searched in nuclear data [Tl], but it is unclear if the 
expectation Q 2 A oc AT A {b) oc A 1 / 3 * 2 / 3 @j is fulfilled at present x. O 2 fT5j. 



2.2. Theoretical developments 

After earlier studies jSJ > a milestone in the theoretical development 
of saturation ideas in QCD was the MV model |17j . This model treats 
classical radiation from color sources moving ultrarelativistically through a 
large nucleus. With a form for the color correlators in the target and taking 
into account the non-abelian gluon interaction, it gives an explicit form for 
the nuclear gluon distribution in the transverse phase space, 

dNf f d 2 

1 oc 



d 2 bd 2 l 



J Sil. e -U-r [! _ e ~Q 2 s r 2 In [(A QCD r)" 2 ]/4j _ ^ 



Later on, gluon radiation (quantum evolution) of color sources was intro- 
duced, leading to a renormalization group-type equation (the so-called Color 
Glass Condensate [IHj)- The rescattering of the projectile in the nucleus is 
described through Wilson lines whose average on target configurations gives 
the S'-matrix. The key point is that the fields (occupation numbers) become 



2 The fluctuation length of the photon into the qq, oc 1/x, becomes at small x larger 
than the hadron size, so its interaction is coherent. 
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large (Fn V , fg oc l/g 2 ) but the coupling is small, so classical arguments and 
perturbative methods are applicable. A hierarchy of coupled equations for 
n-gluon correlators appears [6*1 US]. 



3. The Balitsky-Kovchegov equation 

In the framework outlined in the previous Section, the BK equation 
[""J [7] is the evolution equation in 1/x for the 2-gluon correlator, decoupled 
from the hierarchy in the iV c — ► oo limit (and with correlations neglected 
OH) 3 . As BFKL, it is infrared stable, mixes all twists and a s is fixed. For 
N(xi,X2) = N X1X2 given by a target average of the Wilson lines of a q and 
a q located at transverse positions x\ and X2 respectively, it reads 

dN xlX2 _ _ fd 2 z (x 1 -x 2 ) 2 rAr 1 
8Y ~ Us J 2tt ( Xl -zf{z-x 2 Y [xiZ+ ZX2 X1X2 iV ^ iV ^' 

with Y = In (x initial /x). Taking into account just the linear terms this 
equation reduces to BFKL. In the DLL limit, GLRMQ is recovered .7]. @ is 
the most simple tool we have at our disposal to study the high-density, non- 
linear regime. Its analytical solutions are unknown. Only some theoretical 
estimates exist (221 1213 1^]- So, numerical methods have been developed. 

These numerical methods usually deal with the situation |6| = \x + 
y\/2 3> |r| = |rc — j/| 4 , either in position space [2II1I21] or in momentum space 
[ID 113 [21 EU] (in this latter case we define <f>{k) = J £pi exp (ik ■ r)N(r)). 
Using different techniques two very interesting properties of the solutions 
have been found. First, function h(k) = k 2 V\(f)(k) gets a constant shape 
in In A; at large Y, moving to the right with constant velocity [21 . Second, 
identifying the position of the maximum of h with the saturation scale Q s , 
the solutions show scaling, i.e. (j){k) = <p{k/Q s ) |271 129j (linked to that in 
r discussed in the previous Section, sometimes called geometrical scaling). 
These two features are illustrated in Fig. H The velocity of the solution 
has been computed, Q 2 oc e^ XY , with A ~ 4.1 -=- 4.6 [2TJ 123 US US EH2|- As 
a consequence of these features, the contribution from the low-A: region is 
reduced with increasing Y, thus offering a possibility to avoid the infrared 
problem of BFKL. These features are independent of the initial condition 
(MV [13, BFKL-like [HJ H5], DGLAP-like j27[, Gaussian 28 or GBW 

EH])- Finally, the shape of the solutions above Q s has been examined 
[3"U] . favoring a log-corrected shape [2H1 over a pure power in k [2*2*] . 

3 It was deduced in the color dipole model [Jj, and also (to my knowledge) in an eikonal 
approach 2"2| an< i m BFKL summing fan diagrams in the large N c limit |21) . 

4 The solution with full b dependence has also been studied numerically {25J and turns 
out to be of great importance to compute the behavior of the total cross section at 
very high energies in this approach. 
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Fig. 2. Left: function h versus k for different y = a s Y = 0, I , 2, . . . , 10. Right: the 
same function for y = 4, 6, 8, f and three different initial conditions, versus k/Q s . 
See |301 for further explanations. 

4. Phenomenology at RHIC 

4-1. Multiplicities and pseudorapidity distributions 

For midrapidity at ^/snn = 200 GeV, 0.2 < m/r/GeV < 10 means 
0.002 < x < 0.1, so x at RHIC may well be too high to apply safely sat- 
uration ideas. Assuming that x is small enough, which we will do in the 
following, multiplicities and their evolution with centrality and pseudora- 
pidity can be computed in saturation, usually within a factorized approach 



dN AB 



dydp T d 2 b pj, 



oc \ ( dk T f c 
Pt J 



(y,M/ Qfl (y,PT-k T ) 



(4) 



with Qa,b the saturation momentum of hadrons A and B at some given cen- 
trality, and pt the transverse momentum. It is still unclear to what extent 
this factorized ansatz is exact and which one is the function / which should 
be introduced into this equation |3J; other approaches in the semiclassical 
framework have also been essayed |32| . With some ansatz for /, a sim- 
ple formula is derived US], ^ oc s x / 2 N part e- x ^ lnfQ^/A^^ q2 _ A i/a^ 
A ~ 0.3 (GBW model). This formula successfully reproduces multiplicities, 
their evolution with centrality and pseudorapidity distributions at RHIC. 
See that the deviation from the scaling with the number of participants 
Npart is due just to the log (coming from a factor l/a s (Q 2 s )), at variance 
with what is found in other approaches j^l] • So the question arises whether 
RHIC data can be explained by initial state effects 5 |3bj . 

5 Elliptic flow in saturation models has also been analyzed |51 135|. 
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4-2. Transverse momentum distributions 

The transverse momentum distribution of partons and particles is ex- 
pected to be suppressed by the presence of a medium, the so-called jet 
quenching. This is usually quantified through the ratios 

dN p A dN AA 
o _ dyd 2 p T d 2 b n _ dyd 2 p T d 2 b /-x 
u pA — -Jaj , r^AA — JTj • y~>) 

JVco " dyd I p T d 2 b IycoU dyd 2 p T d 2 b 

Normalized in this way, this ratio goes to 1 at large pt according to the usual 
QCD factorization expectation, with N co u the number of binary nucleon- 
nucleon collisions. Indeed, such suppression has been observed in AuAu 
collisions at RHIC |36j . Contrary to the jet quenching interpretation, it has 
been argued in [37] that such suppression can be explained by initial state, 
saturation effects, so it should also be present in dAu collisions. These colli- 
sions have been studied at RHIC and an enhancement has been found 38 , 
the so-called Cronin effect measured long ago. This leads to the conclusion 
that the depletion in AuAu is due to final state effects. 

Let us examine in more detail the result in |37]. It has been shown 
|39| that the behavior of / above Q s is crucial to get either suppression 
or enhancement in the ratios The naive expectation is that saturation 
effects are important only for k < Q s . Still, the behavior of / in the region 
k > Q s is driven by evolution H31 EH9 • Non-evolved forms lead generally 
to enhancement |391 14U1 141] . On the contrary, in [37] a form which contains 
evolution features [22] has been used. After numerical studies |3U], it has 
become clear J3U] ESI EI] that quantum evolution does not generate 
enhancement but very efficiently erases any that may be present in the 
initial conditions (see Fig. |3J). Then a prediction is that the Cronin effect 
will disappear at higher energies (LHC) or for forward rapidities in pA, 
corresponding to smaller x in the nucleus. Preliminary BRAHMS data |12] 
in dAu collisions at r\ ~ 3 suggest such effect. Nevertheless, other effects like 
running coupling [281 143j may be important for a quantitative comparison. 

Summarizing, the concept of saturation in small x physics has been 
introduced. Some features of the solutions of the non-linear BK equation 
which arises in this context, have been analyzed. The relevance of saturation 
for the initial stage of a nuclear collision has been discussed. Assuming that 
x at RHIC is small enough to apply saturation ideas, the importance of non- 
linear small x evolution for the interpretation of enhancement or suppression 
of the pt distributions measured there has been shown. 
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Fig. 3. Ratios © in pA (upper plot) and AA (lower plot). In each plot, lines from 
top to bottom correspond to rapidities y = a s Y = 0,0.05,0.1,0.2,0.4,0.6,1.0,1.4, 
2.0. See [HOj for further explanations. 
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